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Abstract 

We consider the effects of varying dispersion and nonlinearity on 
the stability of periodic traveling wave solutions of nonlinear PDE of 
KdV-type, including generalized KdV and Benjamin-Ono equations. 
In this investigation, we consider the spectral stability of such solu- 
tions that arise as small perturbations of an equilibrium state. A key 
feature of our analysis is the development of a nonlocal Floquct-likc 
theory that is suitable to analyze the L 2 (M) spectrum of the associ- 
ated linearized operators. Using spectral perturbation theory then, we 
derive a relationship between the power of the nonlinearity and the 
symbol of the fractional dispersive operator that determines the spec- 
tral stability and instability to arbitrary small localized perturbations. 



1 Introduction 



In this paper, we are concerned with the spectral stability and instability 
of periodic traveling wave solutions u(x, t) = u(x — ct) of a class of scalar 
evolution equations of the form 

(1.1) u t + (-A a u + u p+1 ) x = 0, x,teR, 
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where subscripts denote partial differentiation, u = u(x, t) is a real-valued 
function, and where the pseudodifferential operator A = ^—d^, referred to 
as Calderon's operator, is of order one and is defined by its Fourier multiplier 
as Au(k) = \k\u(k); throughout our analysis, the circumflex denotes the 
Fourier transform taken either on R or an appropriate one-dimensional torus, 
depending on the context. By inspection of the Fourier symbol, we see 
that Calderon's operator can alternatively be defined as A = %d x where T-L 
denotes the Hilbert transform being applied either the line or on the torus, 
depending on the context, and in the x variable. Here, we consider a > \ 
and either p 6 N or p = S with m and n being even and odd natural 
numbers, respectively. 

Equations of the form (|l.ip arise naturally in the modeling of unidirec- 
tional propagation of weakly nonlinear dispersive waves of long-wavelength, 
in which case u represents the wave profile or its velocity and the variables 
x and t are proportional to the distance in the direction of propagation and 
the elapsed time, respectively. In this context, the parameter a > char- 
acterizes the linear dispersion about the zero state; in particular, letting 
c(£) denote the phase velocity of plane waves with frequency £, we find the 
dispersion relation c(£) = for equations of the form (jl.ip . It is the goal 
of this paper to derive conditions on the dispersion parameter a and the 
power p of the nonlinearity for which the small amplitude periodic traveling 
wave solutions of (II. ip are stable. 

Arguably the most common and well studied example of an equation of 
form (jl.ip is the generalized KdV equation 

(1.2) u t + u xxx + ( u p +1 ) x = 

corresponding to a = 2. When p = 1, (|1.2p was proposed by Korteweg and 
de Vries [25] in 1895 to model the unidirectional propagation of surface water 
waves of small amplitude and long wavelengths in a channel. When p = 2, 
(jl.2p is corresponds to the modified KdV equation, which arises as a model 
for large amplitude internal waves in a density stratified medium, as well as 
for Fermi-Pasta-Ulam lattices with bistable nonlinearity. In both of these 
cases, the PDE is completely integrable and hence the Cauchy problem can, 
in principle, be completely solved via the inverse scattering transform. 

Another important class of equations of form (jl.ip arises when a = 1, 
in which case we recover the generalized Benhamin-Ono equation 

(1.3) u t + (-Hu x + u p+1 ) x = 0. 

When p = 1, Benjamin [5] and Ono [28J independently derived (jl.3p as a 
model for the unidirectional propagation of internal waves in deep water. 



2 



In this case, (jl .3[) is also completely integrable. Further examples with 
varying values of a can be derived in the context of shallow water theory by 
assuming different order relationships between the small quantities 7 := % 
and /3 := tj, where h denotes the depth of the water at rest, and a and A 
denote characteristic amplitudes and wavelengths of the waves searched for, 
respectively, corresponding to modeling in the small amplitude and small 
wavelength regime. Of particular interest, although outside the scope of 
our analysis, we point out that in the cas^] a = — ^, (jl.ip was recently 
shown by Hur [19] to approximate up to quadratic order the surface water 
wave problem in two spatial dimensions in the infinite depth case, hence 
generalizing Whitham's equation in [32] . 

The stability of the solitary wave solutions of equations of the form (jl.ip 
have a long and rich history, dating back to the novel work of Benjamin in 
[6] in which the stability of such waves in the case a = 2 was established for 
p < 4. The stability analysis for a > 1 was carried out by Bona, Souganidis, 
and Strauss in |7| where the authors extended the seminal work of Grillakis, 
Shatah, and Strauss [16] to equations of form (jl.ip where, most notably, 
the symplectic form in the Hamiltonian structure fails to be invertible. In 
[7], the authors establish the nonlinear orbital stability of the solitary wave 
solutions of (jl.ip when a > 1 for p < 2a and the instability of such solutions 
for p > 2a. 

In contrast to its solitary wave counterpart, the stability theory for T- 
periodic traveling wave solutions of equations of the form (jl.ip has received 
considerably less attention, even in the classical case a = 2, and this theory 
is still far from complete. Within this context, stability results typically fall 
into one of two categories: spectral stability to perturbations in C&(R) or 
L 2 (M) (see [HJ [101 [HI El [T7] ) , and nonlinear orbital stability to perturbations 
in L2 er ([0,nT]) for some n £ N (see [B [21 El U EH US CE21 [20]). The 
majority of the nonlinear stability results restrict to the co-periodic case, 
i.e. perturbations in Lp er ([0, nT\) when n = 1, a clearly very restrictive class 
of perturbations, in which case authors are often able to use adaptations of 
the stability theory in [16] to establish orbital stability. The only examples 
the author is aware of that establishes orbital stability for n > 1 are due 
to Deconinck and Nivala |13j [T3] where the authors consider the KdV and 
mKdV equations, relying heavily in both cases on the complete integrability 
of the governing equations. 

Concerning the spectral stability of such solutions, Deconinck and Bottman 
[llj established the spectral stability in L 2 (M) of periodic traveling wave so- 

1 Notice that the operator ^A" is non-singular for all a > —1. 
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lutions of the classic KdV equation (a = 2 and p = 1), while Deconinck 
and Nivala |13|. IT4"] established such a result for the modified KdV equation 
(a = 2 and p = 2). For more general power law nonlinearities, when a = 2 
Bronski and Johnson [8 J analyzed the spectral stability of such periodic trav- 
eling waves of (|1.2p and derived there a geometric index for the stability and 
instability of such waves to perturbations in Lp er ([0, nT]) with n>l, while 
Haragus and Kapitula [TTj established the spectral stability of such waves of 
sufficiently small amplitude when p < 2 and spectral instability when p > 2; 
as mentioned previously, it is this class of small amplitude periodic waves 
that we will be concerned with here. 

It is important to note that all of the periodic stability analyses described 
above are valid only in the local case a = 2. Such results for a £ (0, 2) 
seem to be very few. Most notably, in [3] Angulo and Natali investigate 
the nonlinear stability of periodic traveling wave solutions to equations of 
the form (jl.ip when subject to perturbations with the same period as the 
underlying wave]!. In particular, they establish the nonlinear stability of 
periodic traveling waves of the Benjamin-Ono equation, corresponding to 
a = 1 and p = 1 in (II. ID to perturbations with the same period as the 
underlying wave. As far as the author is aware, this is the only rigorous 
result concerning the nonlinear stability of periodic waves in equations of 
the form (jl.ip when the dispersive operator is nonlocal. Furthermore, the 
author is not aware of any rigorous results concerning the spectral stability 
of periodic waves in equations of the form (jl.ip when a^2. 

It is the intent of the current paper to investigate the spectral stability 
of periodic traveling wave solutions of (jl.ip of sufficiently small amplitude 
when subject to arbitrary small localized, i.e. integrable, perturbations. Our 
main result, stated in Theorem I3.5j establishes that for a > ^, such small 
amplitude periodic wave trains are spectrally stable if p < ^a-x » wnne they 
are spectrally unstable if p > 2 a-i m ^ * s interesting to note here that for 
a > \ these stability conditions are symmetric!! in p and a, suggesting that 
nonlinearity and dispersion play identical rolls in the spectral stability of the 
small amplitude periodic traveling wave solutions of (jl.ip : see the discussion 
directly following the statement of Theorem 13.51 In particular, such T- 
periodic waves always end up being spectrally stable to small T-periodic 
perturbations so that the spectral instability is necessarily of sideband type, 

2 In fact, they consider more general nonlocal dispersive operators than what are con- 
sidered here. 

3 More precisely, the mapping a M> 2a-i defines a continuous involution on (i,oo). 
Interestingly, the unique fixed point of this involution is a — 2 corresponding to the 
classical gKdV equation (|1.2|l . in which case Theorem 13.51 recovers the result of [17] . 
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occurring in Lp er ([0, nT]) for n > 1, and is hence impossible to detect by 
co-periodic stability analyses. 

While our analysis is similar to that given in [17J in the local case a = 
2, relying on a perturbation arguments of the spectrum from the easily 
analyzed constant state, it is complicated by the facts that (1) the existence 
theory no longer follows by elementary phase plane analysis, and (2) the 
absence of a suitable Floquet theory for nonlocal differential equations with 
periodic coefficients, which is necessary to analyze the essential spectrum 
of the linearized operators obtained from linearizing (11. ID about a given 
periodic traveling wave. In the forthcoming analysis, we resolve (1) by using 
a Lyapunov-Schmidt reduction argument similar to that given in [18] in the 
context of the fifth-order Kawahara equation. For (2), we utilize the inverse 
Bloch- Fourier representation of functions in L 2 (M), which is well known 
in the analysis of Schrodinger operators with periodic potentials [29] and 
has been extensively used in the stability analysis of periodic wave trains 
in dissipative systems (see [H [2TJ [22j [231 EH I2Z] and references therein), 
to show that, even in this nonlocal setting, the essential spectrum of the 
linearized operators acting on L 2 (R) can be continuously parameterized by 
the eigenvalues of a one-parameter family of Bloch operators acting on a 
periodic domain. This spectral characterization extends that introduced by 
Gardner [15] in the local setting, and is valid for considerably more general 
nonlocal operators than what is considered here. 

The outline of this paper is as follows. In Section [21 we prove the ex- 
istence and determine asymptotic expansions of periodic traveling wave so- 
lutions of ([Lip of sufficiently small amplitude. The existence argument is 
based on an appropriate Lyapunov-Schmidt reduction argument, the details 
of which are included in Appendix [Aj Section [3] contains our main stability 
results, beginning with a careful characterization of the essential spectrum 
of the linearized operators acting on L 2 (R) in terms of the eigenvalues of a 
one-parameter family of Bloch operators acting on periodic functions. We 
then use spectral perturbation arguments to analyze the spectrum of the 
small amplitude periodic wave by considering the associated Bloch opera- 
tors as small perturbations of those with constant coefficients obtained from 
linearizing about the nearby constant state. As such, the restriction to pe- 
riodic waves of sufficiently small amplitude is essential in our argument. In 
particular, our analysis gives no information about the stability or instabil- 
ity of periodic waves with large amplitude when subject to small localized 
perturbations. Finally, we conclude with an appendix in which we give the 
proofs of the existence result and the Bloch-wave decomposition. 
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2 Existence of Periodic Traveling Waves 



In this section we analyze the set of periodic traveling wave solutions of (jl.lj) 
of the form 

u(x, t) = u{x — ct), c£l, 

where the function u(-) is a real-valued periodic function of its argument. 
Due to the scaling properties of (jl.ip . we can with out loss of generality 
assume that c = 1 in which case such solutions of (II, ip arise as stationary 
solutions of the PDE 

(2.1) u t + (-A a u-u + u p+1 ) x = 0, 

or, equivalently, after performing a single integration, solutions of the non- 
local profile equation 

(2.2) - A a u - u + u p+1 = b 

where b £ R is a constant of integration taken to be in general non-zero. 

We begin by considering the equilibrium solutions of (|2.2p for \b\ -C 
1. Notice when 6 = 0, there are in general two non-negative equilibrium 
solutions u = and u = 1. In the classical case when a = 2, it follows 
by straight forward phase plane analysis that in the (u, u') phase plane, 
u = is a saddle point associated with a homoclinic orbit (solitary wave), 
while the equilibrium u = 1 is a nonlinear center. Thus, when a = 2 and 
6 = it is clear that there exists a one-parameter family of periodic orbits 
{(u 7 , ^)} 7e r j\ °f period T 7 such that 

lim T 7 = , lim T 7 = +oo. 

Moreover, these qualitative features persist for \b\ <^ 1. 

When a 7^ 2, however, this classical picture breaks down as we can 
no longer make sense of the phase space in the same way. Nevertheless, for 
\b\ <g 1 there exists a unique equilibrium solution u = Qt of (j2.2j) continuing 
from Qo = 1; indeed, a straightforward calculation provides the expansion 

Q fe = i + -6-^& 2 + o(|6| 3 ) 

p 2p A 

valid for |6| <C 1. Seeking nearby periodic solutions of (|2.2p . we set u(x) = 
P(kx), where P is 2-7r-periodic and k > denotes the wave number, and 
require that P and k satisfy the rescaled nonlocal profile equation 

-k a A a P - P + = 6 
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posed on a 27r-periodic domain. Here, we consider for each s > the operator 
A s as being defined on the torus. In particular, we consider A s as a closed 
operator on L 2 (M/2vrZ) with dense domain H s (R/2ttZ) being defined via 
Fourier series as 

A s /(x)= Yl \ k \ SelkX f(k), s>0. 
fcez\{o} 

It is clear from this definition that A s is invertible for each s > when 
restricted to the mean-zero subspace of H s (M./2ttX), and we define its inverse 
A~ s on this subspace via Fourier series as 

A-7(*)= E 

k<=Z\{0} 1 1 

see |30j for a recent interesting discussion of A s acting on the torus. With 
these definitions, for a > \ one can use a Lyapunov-Schmidt reduction 
to verify the existence of a two-parameter family of small amplitude even 
periodic traveling wave solutions u a b of (|2.2|) existing in a neighborhood of 
the equilibrium solution u = Qb\ see Appendix [A] for details. 

In the parameterization of the periodic solutions u a> b(x) = P a ,b{k a ,bx) of 
(|2.2p given by Theorem lA.il the functions P a ^ are 2-7r-periodic even solutions 
of the rescaled nonlocal profile equation 

(2.3) - K )b A a v -v + v p+1 = b 

such that Po,fe = Qb and k^ b = (p+l)Q^— 1; here, we have used the fact that 
k a ,b > for | (a, b)\ <C 1. Furthermore, the parameter a is precisely the first 
Fourier coefficient P a ,b-> and we have the relations Pa^z + ir) = P- at b(z) and 
k a fi = k_ a b valid for all |(a, b)\ <C 1. Taking into account the translation 
invariance of (12. 2p . it follows that for a fixed a > ^ we can find a three- 
parameter family of small amplitude periodic traveling wave solutions of 

m- 

A key feature of the functions P a ^ and k a fi described in Theorem lA.ll is 
that these depend analytically on a and b for |(a, 6)| <C 1. The next lemma 
exploits this fact to provide us with explicit expansions of the functions P a> b 
and k a b valid for sufficiently small a and b. These expansions are crucial in 
the forthcoming stability analysis. 



7 



Lemma 2.1. For sufficiently small o,f) 6 R and a > \, the functions P a ^ 
and k a £ in Theorem \A.1\ can be expanded as 

P a ,b(x) = Q b + cos(z)a + Q cos(2z) - 1^ a 2 + 0(\a\{a 2 + 6 2 )) 

Proof. Taking 6 = for now and recalling Theorem IA.ll we begin with the 
following small amplitude (Stokes wave) ansatz for P a $ and k Qj o 

Pa,o( z ) = 1 + acos(z) + a 2 v 2 {z) + a 3 v 3 (z) + 0(a 4 ) 
k%=p + a 2 k 1 + 0{a A ) 

where each Vj is even and 27r-periodic in the z variable. Substituting these 
expansions into the profile equation (|2.3p yields a hierarchy of compatibility 
conditions. The 0(1) equation is trivially satisfied, while the 0(|a|) equation 
reads 

(-A a + l)cos(z) = 
which again holds. The 0(a 2 ) equation now reads 

-A a v 2 + v 2 = - gfc-ti] (i + cos (2z)) 

which has for any constant A £ R the even solution 

V2 ( z ) = Pk±H Qcos(2z) - + Acos(z). 

Notice, however, that we must take A = by the definition of the parameter 
a. Continuing to the 0(a 3 ) equation, we find 

p (-A" + l)v 3 = h cos(z) - p ( p +^ 2 (cos(z) + cos(3z)) + p{p + ^ cos ( z ) 

p{p- l)(p + 1) 
27 ( 3 cos (^) + cos(3z)) , 

which is readily seen to have a resonant solution containing a term propor- 
tional to zsin(z) unless 

_ p(;p+l) 2 p(p+l) 2 3p(p-l)(p+l) _ 
1 24 4 24 
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i.e. unless k\ is chosen so that 



p(p + l)(p + 4) 
1 12 

This completes the expansions to the desired order when 6 = 0. The expan- 
sions when \b\ -C 1 are obtained analogously. □ 



3 Spectral Stability to Localized Perturbations 

We are now ready to begin discussing the stability analysis of the small 
amplitude periodic solutions P a ,b{ka,b ) constructed in the previous section 
under the flow induced by the PDE (|2. 1|) . As mentioned in the introduction, 
we are interested here in the spectral stability of such waves when subject 
to small localized, i.e. integrable, perturbations on R. In the classical case 
when a £ 2N, the associated spectral problem obtained from linearizing 
about a given wave P a ^ is that for an ordinary differential equation with pe- 
riodic coefficients. As such, it is an easy calculation to see that the spectrum 
of this linearized operator, considered as an operator on L 2 (R) is purely es- 
sential and agrees with the continuous spectrum. In this classical 
common method for analyzing the essential spectrum of the linearization 
is to use Floquet theory to provide a continuous parameterization of the 
essential spectrum by the eigenvalues of an associated one-parameter fam- 
ily of linear operators considered with periodic boundary conditions; see 
P[lQl[T7], for example. 

When a > ^ is not an even natural number, however, the linearized 
operator associated with P a ^ is nonlocal and hence the classical Floquet 
theory does not apply. Nevertheless, we find in the next section that we 
can still reduce the spectral stability problem on L 2 (M) to a one-parameter 
family of eigenvalue problems with periodic boundary conditions. This char- 
acterization of the essential spectrum is given by utilizing the Floquet-Bloch 
transform defined on L 2 (R). Once this characterization is established, we 
use an appropriate spectral perturbation theory to analyze the eigenvalues 
of the associated one-parameter family of Bloch operators when |(a, 6)| <C 1 



3.1 Characterization of the Essential Spectrum 

To analyze the spectral stability of the small amplitude periodic waves ob- 
tained in the previous section, we fix a > ^ and set z = k a f,x and s = k a ^t 
in flUQ to get 

v s + (-\k a , b \ a A a v-v + v p+1 ) z , 
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where now A acts on R in the z-variable. Linearizing about P a ,b and consid- 
ering solutions of the linearized equation of the form v(z,t) = e xt v(z), with 
A G C and v(-) G L 2 (R), leads to the spectral problem 

M a fiV := d z C afi v = Xv 

considered on H a+1 (R), where here A denotes the spectral parameter and 

C a ,b :=\k a:b \ a A a + l-(p+l)Pl b 

is considered as a closed, densely defined operator acting on L 2 (R). In this 
case, spectral stability is defined by the condition that the operator Ai a ,b 
have no spectrum in the open right half plane. However, however, that 
the symmetries of the spectral problem imply the spectrum is symmetric 
with respect to the real and imaginary axes, and hence spectral stability is 
equivalent with all the spectrum of the operator Ai a ^ being confined to the 
imaginary axis; we will discuss this more in the next section. 

In order to analyze the spectrum of the linear operator M a ,b acting on 
L 2 (R), we recall now some facts about the Floquet-Bloch decomposition of 
L 2 (K); as mentioned above, the standard Floquet theory does not apply for 
a $l 2 N since the spectral problem for M a ,b is n °t in the form of an ordinary 
differential equation. To this end, notice that given any v G L 2 (R) we can 
express v in terms of its inverse Bloch representation as 

v(x)= I''* e* x v{Z,x)d£ 

J -1/2 

where v(£, x) := Ylkez e tkx v(£ + k) are 27r-periodic functions of x and where 
v(z) := J R e~ lu)Z v(z)dz denotes the standard Fourier transform of v. In- 
deed, the above formulas may be easily checked on the Schwartz class by 
grouping frequencies which differ by one in the standard Fourier transform 
representation of v: 

r l/2 fl/2 

v(z) = V / e i{k+ ^ z v{k + j)dk = / e ikz v(£,z)dz. 
jeZ J-i/2 J -1/2 

The Bloch transform B : L 2 (R) L 2 ([-l/2, 1/2); L 2 (R/2ttZ)) given by 
B(v)(^,x) := v(^,x) is then well defined, bijective, and continuous; in fact, 
using the classical Parseval theorem we find for all v G L 2 (R) that 

rl/2 f-2-rr 

(3.1) \\v\\ 2 L2{R) = 2tt / / \B(v)(Z,z)\ 2 dzd£ 

J — 1/2 J 
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so that the rescaled Bloch transform \^2ttB is an isometry on L 2 
Furthermore, we find given any v G L 2 (R) that 



B(M a>b v) (t,x)=M aM m,-))(x) 

where M a M ■ H a+1 (R/2irZ) C L 2 (M/2vrZ) -»• L 2 (M/2vrZ) is a defined by 

M aM := e-* x M a ,be* x , £ G [-1/2, 1/2). 

The operators M. a ,b£ are called the Bloch operators associated to M a ,b> 
and from above they correspond to operator-valued symbols of A4 a ,b under 
B acting on L 2 (IR/2-7rZ). The next result relates the spectrum of M. a ,b 
acting on L 2 (R) to that of the corresponding one-parameter family of Bloch 
operators {•M.a,b,^}^\-U2 1/2)' providing us with a Floquet-like theory that 
is suitable for our needs. 

Proposition 3.1. Consider the operator Ai a j, acting on L 2 (R) with domain 
H a+1 (R) and the associated Bloch operators {■A / la,6,f}^ e r_i/2 1/2) acting on 
L 2 (R/2ttZ) with domain H a+1 (R/27rZ). Then for any A G C, the following 
statements are equivalent: 

(i) A belongs to the spectrum of the closed operator M a b acting on L 2 (IR). 

(ii) There exists a £ G [—1/2,1/2) such that A belongs to the spectrum of 
the closed operator -M a ,6,£ acting on L 2 (M/2-7rZ). 

(Hi) There exists a nonzero function V G L 2 (R/2tt7j) of the form V(z) = 
e^ z v{z) for some £ G [-1/2,1/2) and v G H a+1 (R/2TrZ) such that 
(M a , b - AI) V = 0. 

The proof of Proposition 13.11 is contained in Appendix |Bj Using this 
result, it follows that 

°~L 2 («) (M a ,b) = |J Crz,2 ef ([o j27 r]) (M a ,b,l-) 

?e[-i/2,i/2) 

so that, in particular, the Bloch transform provides a continuous parametriza- 
tion of the essential spectrum of the operator M. a ,b by the discrete spectrum 
of a one-parameter family of Bloch operators. As a result, rather than 
analyzing the essential spectrum of the operator A4 a ,b directly, we can in- 
stead choose to study the point spectrum of the operators M a ,b,£ for each 
£ G [—1/2, 1/2). This investigation is the subject of the next section. 

4 This is in fact a special case of a more general class of generalized Hausdorff- Young 
type inequalities, following from interpolating (|3.1[) with the triangle inequality, that are 
known to be satisfied by the Bloch transform; see [24] for details. 
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3.2 Analysis of the Unperturbed Operators 

We begin by considering the stability of the constant state Po,o = Qo = 
1. Later, we will treat the linearized operators about the nearby periodic 
solutions as small perturbations of those linearized about Qo. Indeed, it is 
straightforward to establish the estimate 

\\M aM -M , ,d = O(\a\ + \b\) 

as (a, b) — > (0,0) uniformly in the Bloch parameter £ G [—1/2,1/2). A 
standard perturbation argument then guarantees the spectrum of M a ,b,^ 
and .Mo,o,£ stay close for sufficiently small (a,b). More precisely, we have 
the following result. 

Lemma 3.2. Let p > 1 and a > ^. For any 5 > i/iere exists an e > 
stic/i i/iai, /or any £ G [—1/2, 1/2) and any (a, 6) £ M 2 mt/i ||(a, 6)|| < e, the 
spectrum of Ai a ,b£ satisfies 

a (Maw) C {A G C : dist (A, a (M 0M )) < 5} . 

We now analyze the spectrum of -Mo,o,§ posed on L 2 (M./2irZ). Since 
.Mo,o,£ has constant coefficients, we find by a straightforward Fourier analysis 
argument that 

<r(M 0t o£) = {X = iu n £ : n G Z} C iR 

for each fixed £ G [—1/2, 1/2), where the uj n £ are determined by the linear 
dispersion relation 

(3.2) oj(k) := kp (\k\ a - 1) 

through cj n £ := oj(n + £). Notice that every A G a(Mo,o,^) is a semi-simple 
eigenvalue with algebraic and geometric multiplicity given by the number of 
distinct n£Z such that A = ito n £ with associated eigenfunction e n := e mz . 

To study the behavior of these eigenvalues for small (a, b) G M 2 , notice 
that the spectrum of the operator M a ,b is symmetric with respect to both the 
real and imaginary axis. Indeed, since the coefficients of M a ,b are rea l valued 
it follows that its spectrum is symmetric with respect to the real axis. In 
terms of the Bloch operators this implies that a{M. a ^^) = a{M. a b _^). Fur- 
thermore, noting that M a ,b anti-commutes with the isometry S : L 2 (R) — >• 
L 2 (R) given by 

Sv(z) = v(—z) 

we find that the spectrum of A4 a ,b is symmetric with respect to the origin. 
It follows that the Bloch operators satisfy M. a ,b^S = —SM. a ^^ so that 
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&(M a ,b,£) = —&(M a ,b-s)- Finally, recalling that k a ^ is even in a and that 
P a ,b{z + vr) = P_ atb (z) for all |(a, b)\ «lwe see that a(M a ,b) = v(M- a ,b) 
and that <r{M a ,b,{) = cr(M- a ,b^)- As a result, we find that the spectrum of a 
given Bloch operator M a fi£ is symmetric with respect to the imaginary axis. 
It follows then that when eigenvalues of M a ,b,^ bifurcate from the imaginary 
axis they must bifurcate in pairs resulting from collisions of eigenvalues on 
the imaginary axis. A well known result from the study of Hamiltonian 
systems tells us that when two purely imaginary eigenvalues collide, the 
collision will not result in a pair of eigenvalues bifurcating from the imaginary 
axis provided both eigenvalues have the same Krein signature; see |33j . 

We now consider the location of the eigenvalues more carefully, in par- 
ticular watching for sets of eigenvalues that collide for a fixed £. Notice 
by the symmetry property a (Ai a ,b,$) = a {M-a,b-i ) we ma Y restrict our 
consideration to Bloch frequencies £ G [0,1/2]. Now, when £ = we find 
that 

w_i,o = w o,o — ^i,o = 

and 

. . . < W_ 3j o < < < w 2>0 < uj 3fi < .... 

Furthermore, for £ G [0,1 /2] we have 

U [2p(2 a - l),oo) , |n|>2 
and that 

■^ c [I(f- 1 ) , t (GO" -1 )]' N£1 ' 

which naturally provides us with a spectral decomposition 
° (Mom) = a ± (Mom) u °"2 (M ,o,t) 

for -A/fo,o,§ with 

fi (-^0,0,^) = {A G C : A = ou n £ for some \n\ > 2} 
0-2 (M ,o,t) = {^-1,5,^0,5,^1,5} , 

with the property that for any v in the infinite dimensional spectral subspace 
associated with a% (A^o,0,$) satisfies 

(C 0M v,v) > (Y|) -l)p\\v\\ 2 
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uniformly for £ G [0, 1/2], where 




is considered on L 2 (IR/27rZ); in particular, it follows that all eigenvalues in 
ci(A^o,o,^) have positive Krein signature for £ G [0, 1/2]. By a standard per- 
turbation argument, we find that the above properties persist for sufficiently 
small a and b. More precisely, one has that for a and b sufficiently small we 
have a spectral decomposition 



for | (a, 6) | <C 1 and £ G [0,1/2], where the spectral subspace associated 
with o"2 (At 6 is three dimensional, and the infinitely many eigenvalues 
associated with o\ [M. a ,b,i) all have positive Krein signature. Notice this 
latter property implies that all the eigenvalues in o\ (M. a ,b,e) are purely 
imaginary for |(a, b)\ sufficiently small. 

It now remains to determine the location of the eigenvalues in a 2 {M a ,b,^)- 
These are smooth continuations for small a, b G M. of the eigenvalues iuj-±£, 
iojQ^, and of M a ,b,£- First, notice that two such eigenvalues collide if 
and only if £ = 0, when ^±1,0 = wo,o = 0. Furthermore, for any £0 £ (0, 1/2) 
there exists a constant cq > such that 

- u k ^\ > c p, j, k G {—1, 0, 1}, j^k, £G [£o,l/2]. 

Consequently, these eigenvalues will remain simple and distinct under small 
perturbations (for sufficiently small a and b) for all £ G [0, 1/2]. In particular, 
it follows that for any £0 G (0, 1/2] we have 



for sufficiently small a and b. Thus, we have reduced the problem to locating 
the eigenvalues iw±\£ and zujq^ for a, b, and £ small. 

3.3 Location of the Critical Bloch Spectrum 

To this end, for |(a, b, £)| < 1 we compute a suitable basis for the spectral 
subspace associated to (T2{M. a ,b,t) and compute the 3x3 matrix 



o {M aM ) = a x {M aM ) U a 2 {M aM ) 



such that 



°\ (M aM ) n a 2 (M aM ) = 0, 



1 (M 



a. 



, M )ciR, £G[£ ,l/2] 




(Vj{a,b, £),% (a,6,£)) 



,M a ,b,tVk(a,b,0 



> 



) 



-I j,fc=l,2,3 
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representing the action of M a ,b,£ on this space. First, at a = b = we 
have that the operator .Moo£ has constant coefficients and a basis for the 
critical spectral subspace in this case is spanned by the functions 1 and 
e ±x associated with the eigenvalues iojq^ and ioJ±x£, respectively. Working 
instead with the real basis 

7/0 = cos(z), 771 (z) = sin(z), r] 2 (z) = 1 

a direct calculation yields 

/ § + w_i >f ) \ 

\ iwo,( / 

valid for any £ G [0, 1/2]. Here, we are using (/, g) = J* Q 27r f(z)g(z)dz. 

Next, at £ = we claim that A = is an eigenvalue of M a ,b with 
algebraic multiplicity three and geometric multiplicity two. Indeed, notice 
that since (II. lj) is invariant with respect to spatial translations it follows 
that 

M a ,b,od z P a>b (z) = 

so that A = is indeed an eigenvalue of M a ,b,o- Furthermore, differentiating 
the profile equation (|2.3p with respect to the parameters a and b yield 

M a , b , d a p a , b (z) = -d a (\k a , b \ a )A a d z p a , b 
M aA0 d b p ajb (z) = -d b {\k a ^ b \ a )K a d z p^ b 

so that, in particular, we find 

M a , b ,o (d b (\k a , b \ a ) d a P a , b (z) - d a (\k a>b \ a ) d b P a>b (z)) = 0, 

giving a second function in the kernel of M a ,b,o- Finally, by a straight 
forward computation, using (12.31) and the fact that J~- a ,b,od z P a ,b = 0, we find 

M aA0 Pa,b( z ) = ~P (IM aAa + !) d * p a,b 

so that 



M a , bfi (d b (\k a!b \ a )P a>b -p\k a>b \ a d b P aib ) = -pd b (\k a , b \ a )d z P a>b , 

giving a 2-7r-periodic function in the generalized kernel of M a ,b,o ascending 
above the translation mode d z P a ^. In particular, this shows that for all 
values of a and b sufficiently small, the three eigenvalues iloq^ and iu±i£ 
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all vanish when £ = 0, and that the eigenvalue A = indeed corresponds to 
an algebraically triple and geometrically double eigenvalue of the operator 
M a ,b,o acting on L 2 (M/2vrZ). 

Using the expansions in Lemma [2. II we now obtain a basis for the critical 
eigenspace which is compatible with the basis found at a = b = 0. In 
particular, we take 

% = "4t CM") d a Pa,b( Z ) ~ da CM") 9 b P a ,b^)) 
p+l 

2 p — 1 p+l 2 

= cos(z) a H cos(2z)a cos(2;)6 + 0(a 2 + b 2 ) 

6 6 p 

m = --d z P ab = sm{z) + ^^sm(2z)a + 0(a 2 + b 2 ) 
a b 

V2 = d b {\k a)b \ a ) P a<b - p\k a)b \ a d b P afi 

= 1 + (p + 1) cos(z)a - ^-^6 + C(a 2 + b 2 ). 

P 

In this basis, a straightforward calculation yields 


B a ,b,o = | CJ a>b 




where 



0"a,6 = : LJ r = £>a<9 fe ( fca,6 J 



= p(p + l)a - 2(p + l)a6 + C(|a| (a 2 + 6 2 ) 

Furthermore, in this basis we note that the projection of the identity matrix 
can be expanded as 



(rj j (a,b,0),7] k (a,b,0)) 
(r] j (a,b,0),r] j (a,b,0)} 



1 2±^a 

1 | + 0{a 2 + b 2 
j,fc=i,2,3 \ 1 



In order to analyze the breaking of the algebraically triple eigenvalue at 
A = of the operator M a ,b,e, for |£| <C 1, we need the following result which 
allows us to expand the nonlocal dispersive operator d x A a with respect to 
the Bloch frequency £. 
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Lemma 3.3. Let a > and f G H a+1 (R/27rZ) be fixed. Then for all 
|£| <C 1, we /ioue 

e-^d x A a e^f(x) = d x A a f(x) + %\£\ a f{0) 



+ (E i ( 2 a / t 1 1 ) Aa " 2( ^ 1) ^ 1 ) v K x ) 



+ E 



v r=l 



a + 1 
2r 



where (™) := m ( m 1 )( m ^"^ m n+1 - > denotes the generalized Binomial coeffi- 
cient, defined for m G C and n G N, and P denotes the orthogonal projection 
of H a+1 (R/2irZ) onto the subspace of mean-zero functions. 

Proof. Using the Fourier series representation of / we find 

e-* x d x A a e* x f(x) = ^i(k + £)|fc + £\ a e ikx f(k). 
kei 

Since |£| <C 1, we have |/c + £| = |fc| + sgn(/c)£ so that, for k ^ 0, the term 
| A; + £\ a may be expanded using Newton's Binomial series as 

i fc +^r = E ( )\k\ a - m sm(krr, 

and hence 

oo 

(k + 0\k + CI" = |fc| a * + E sgnCfc)" 1 " 1 



m=l 



n 



77? 



+ 



m — 1 



|A;| 



a— (m— 1) £ m 



valid for 7^ 0. If m = 2Z — 1 for some / G N, then 

sgn(A:) m " 1 |A:| Q - (m - 1) = \k\ a ~ 2 ^ 

while if m = 2r for some r G N, we have 

sgn(A;) m - 1 |A ; | a -( m - 1 ) = k\k\ a ~ 2r . 

Using Pascal's rule (™) + ( n "\) = ( m + 1 ), valid for all m G C and n G N, we 
find for k 7^ 



00 / -1-1 

(fc+oi*+fr = *ifci a +E(2z_i 



+E 

r=l 



a + 1 
2r 



| fc |a-2(i-l)^2i-l 
k\k\ a - 2r C 2r , 



from which the claimed expansion follows. 



□ 
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Remark 3.4. Notice that, for a given a > the smoothness of the operator 
e~ i * x d x A a e i € x in £ near £ = is determined by the smoothness of the term 
near the origin. In particular, for a 6 2N one obtains an analytic 
expansion, while for positive a £ 2N one obtains only a cL a J +1 expansion. 
It is expected then that for every a > the spectrum of M. a ,b,£ bifurcating 
from the (A, £) = (0, 0) will be at least C 1 in a neighborhood of £ = 0. 

From Lemma I37LI we can expand for |(o, b, £)| <C 1 the operator A4 0) 6^ 
to any desired order in £. For our purposes, it is sufficient to identify the 
term of first order in £, which is 

= i((a + l)k« b A a + l-(p+l)Pl b ). 

Notice that outside the small amplitude regime |(a, b)\ <C 1 one would have 
to consider higher order terms in the expansion as well. Indeed, in the 
classical case a = 2, corresponding to the KdV equation, one typically must 
expand the Bloch operator to second order in £ ; see [10] . For large amplitude 
periodic waves then, Remark 13.41 then seems to indicate one must take care 
to consider the cases a > 1 and a £ (0, 1) separately. In our case, however, 
we will see below that the restriction to small amplitude solutions requires 
only an expansion to first order, hence bypassing this added difficulty. 

Continuing, we note that the basis {r]o, rji, 772} defined above can be 
extended to a basis of the three-dimensional eigenspace bifurcating from the 
generalized kernel of A4 a ,b,o f° r sufficiently small a, b and £. This provides 
us with an expansion of the form 
(3.4) 

B aM = Boflt+B^o + aZBx + btlh + O (|£|(a 2 + b 2 ) + |e| min ( 2 ' Q+1 )(a + 6)) . 

In order to determine the three eigenvalues of B at b£, we consider the char- 
acteristic polynomial 

D(\) = - det (B aM - XI) = A 3 + c 2 A 2 + c x A + c 

where the coefficient functions Cj = Cj(a,6,£), defined for |(a, 6, £)| <C 1, 
depend smoothly on the parameters a, b and are C 1 in £. Analyzing the 
dependence of the Cj on a, 6, and £ more carefully, notice that since the 
spectrum of M a b{, is symmetric with respect to the imaginary axis, consid- 
ering the coefficients of D as symmetric functions of its roots, we see that the 
C2 and Co must be purely imaginary whereas c\ must be real. Furthermore, 
since 

O (M a ,b,t) = O- (M-aM) » a ( M aM) = ° (M a ,b,-t), 



(3.3) Tr c M aM 
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it follows that the coefficients Cj must all be even in a and that c 2 and cq are 
odd in £ while c\ is even in £. Together with the expansion of £> a ,6,? & bove, 
these properties imply that the polynomial D can be expressed as 

DW = -(«)»« (^a^- A) 

so that the critical eigenvalues of .Ma^f bifurcating from the (A,£) = (0,0) 
state are of the form A = ip£X where X is a root of the cubic polynomial 

Q(X) = - det (J^ B aM -X^j =X 3 + d 2 X 2 + d x X + d , 

where the coefficients dj are real- valued and even in a and £. To determine 
whether the roots of Q, and hence the three critical eigenvalues, lie on the 
imaginary axis or not we consider the discriminant 

A aM = 18d 2 did + d\d\ - Ad\d Q - Ad\ - 27d%, 

which here, by the symmetry properties of the coefficients dj, can be ex- 
panded near for |(a, b, £)| <C 1 as 

(3.5) A aM = A 0M + 7 a 2 + ((a 2 (a 2 + \b\ + £ 5 )) 

for some appropriate 5 > 0, chosen independently of £, and some constant 
7 G R still to be determined. 

Notice that the operator A^o,fe,g has constant coefficients and is clearly 
skew-adjoint. Consequently, the spectrum of Mo,b,^ must lie on iM. so that, 
in particular, Aq 55 > 0. More precisely, using (j3.2j) these eigenvalues are 
given explicitly by 

iz\k , b \ a (ier - 1) , i(±i + oim q (i ± 1 + er - 1) . 

It follows that 

= L Sr (Mb <x)M-Z; a)A 2 {t- a)) 2 , 

where, for |£| <C 1, we have the expansions 

A X & a) = -1 + (1 + er +i - m a = (« + m + o(ici 1+£ ) 

4 2 (£; a) = -2 + (1 - C) Q+1 + (1 + £) Q+1 = a(a + l)f + 0(\i\ 2+£ ) 
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for some appropriate e > chosen independently of £. Thus, for |£| <C 1 we 
have 

A 0,b,5- ~6 + <-Al£l ), 

which is clearly positive for £ sufficiently small. 

To compute the coefficient 7 in (|3.5|) . notice it is enough to compute 
Aa,o,o to 0(a 2 ). To this end, using (|3.3p we find that the matrix £>i in (|3.4p 
is given by 



01 



/ (a-l)p(p + l) + ^zl) 



(q-l)p(p+l) | p(2p-l) q q 



which gives, using Mathematica [2B], the expansion 

. a(a + 4)(1 + a) 4 . 1S /\ 2a — l . ., , 
A«,o,o = — ^ ~(p + l) I - -r; — P a + °( a )■ 



4 + a 

Together, the above analysis establishes our main result. 

Theorem 3.5. Let a > \ and p > l be fixed such that either p £ N or 
p = — where m and n are even and odd integers, respectively. Then the 
small amplitude periodic traveling wave solutions u a ^ = P a ,b(k a ,b') of (jl.ip 
constructed in Theorem \A.l\ for \(a,b)\ <C l are spectrally stable if p < 2a -1 > 
and are spectrally unstable if p > 2a— 1 • 

Remark 3.6. As mentioned in the introduction, the spectral instability de- 
tected in Theorem I3.5I is of sideband type, and is hence not detectable if one 
restricts to perturbations with the same period as the underlying wave. 

As far as the author is aware, Theorem I3.5I is the first rigorous spec- 
tral stability/instability result in the periodic context for nonlocal KdV like 
equations of form 111. IB . It is interesting to note that the conditions between 
p and a in Theorem I3.5I to establish stability or instability are symmetric. 
More precisely, the mapping a h-> 2a-\ defines a continuous involution on 
(|,oo), so that for a > | we have 

4 + a 4 + p 

p < <==^ a < 



2a -I 2p-l 

Thus, in the case of small amplitude periodic traveling wave solutions of 
(jl.ip . dispersion and nonlinearity seem to play identical rolls, at least as far 
as the stability analysis is concerned. 
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Notice that when a = 2, corresponding to the classic gKdV equation 
(jl.2p . Theorem 13.51 recovers the result of Haragus and Kapitula in |17| . 
yielding spectral stability for p < 2 and spectral instability for p > 2. In- 
terestingly, this classical case corresponds to the unique fixed point of the 
involution discussed in the preceding paragraph. As one begins to decrease 
a from this classical case, the critical power of the nonlinearity which signals 
a switch in stability increases, and tends to infinity as a — > | from above. 
This behavior is exactly the opposite of what occurs in the solitary wave 
case, where the critical power decreases linearly with the dispersion param- 
eter a as a is decreased; see [Tj- This illustrates an interesting difference 
between the solitary and periodic theories. In particular, notice that the 
restriction here that p > 1 implies that all such small amplitude periodic 
waves are spectrally unstable to localized perturbations when a > 5. 

Another model of particular interest in our analysis is the generalized 
Benjamin-Ono equation (|1.3p . corresponding to a = 1. In this case, Theo- 
rem [33] implies that for \b\ <C 1, the small amplitude periodic traveling wave 
solutions of (|1.3p are spectrally stable for p < 5 while they are spectrally 
unstable for p > 5. 

It is also interesting to analyze Theorem 13.51 by fixing the nonlinearity p 
and asking what dispersion parameters a leads to spectrally stable waves, 
although these are clearly related to the results for fixed a by our previous 
comments. For p = 1, we find that such small amplitude periodic traveling 
waves of (jl.ip are spectrally stable provided a < 5, while they are spectrally 
unstable if a > 5. For p = 2, corresponding to a dispersion generalized 
modified KdV equation, Theorem 13.51 gives stability for the small amplitude 
solutions u a fy for a < 2 and spectral instability if a > 2. Furthermore, for 
p = 3 we find such solutions are spectrally stable for a < g and for p = 4 
they are stable if a < |. In particular, we find that for a fixed p > 2, one 
can find small amplitude periodic traveling wave solutions of (jl.ip that are 
spectrally stable to small localized perturbations provided the dispersion 
strength is chosen to satisfy | < a < < 2. This seems to be an 

interesting result complimenting the spectral instability results of Haragus 
and Kapitula [T7j in the local a = 2 case. 
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Appendix 

A Existence of Small Periodic Traveling Waves 

In this appendix, we utilize a Lyapunov- Schmidt reduction to establish the 
existence and basic regularity properties of small amplitude periodic trav- 
eling wave solutions of the nonlocal profile equation (|2.2p . The method of 
proof parallels that given for Theorem 1 in [18] in the context of the fifth 
order Kawahara equation. 

Theorem A.l. Let a > \, andp > 1 be fixed such that either p G N or else 
P = ~ w ^ ere 171 and n are even and odd integers, respectively. Then there 
exists constants ao,bo G R such that for any fixed b G (— bo,b ) the nonlo- 
cal profile equation (|2.2p admits a one-parameter family of even, periodic 
solutions W,fc}ae(-ao,a ) of the form 

U a ,b(x) = Pa,b(ka,bx) 

where P a< b is 2ir-periodic and smooth in its argument. Moreover, the follow- 
ing properties hold: 

(i) The map k : (— ao,ao) x (—bo, bo) — > R is analytic, even in a, and 
satisfies 

k« b = k*(b) + ~k(a,b) 
where k*(b) = (p + — 1 and 

k(a,b) = y^k2n(b)a 

n>l 

for any \(a,b)\ <C 1 and some positive constants Kq and po > ao- 

(ii) The map (— ao, ao) x (— bo, bo) B (a,b) h-> P a ^ G H a (JH/2ir'L) is analytic 
and can be expanded as 

P a ,b(z) = Q b + acos(z) + Y, Vn, m (b)e l(n - m)z a n + m 

n,m^0, n+m>2 
n— m^±l 

where p n , m G R are such that p n ,m(b) = Pm,n(b) with 

Co 

\Pn,m(b)\ < 

Po 

for any \b\ <C 1 and some Co > 0. 



k 2 n(b) 



< 



AO 
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(Hi) The Fourier coefficients p n (a,b) of the 2ir -periodic function P a ,b> 

Pa,b(z) = Y J Pn(a,b)e mZ 

are real and satisfy po(a,0) = 1 + 0{a 2 ) and p n (a,0) = 0(|a| n ) for 
n ^ as a — > 0. Moreover, the map a \— > p n (a,b) is even (resp. odd) 
for even (resp. odd) values of n. 

Proof. Renormalizing the period to 2n, let k G R + denote a wavenumber 
and set z = kx so that the rescaled profile equation becomes 

(A-l) - k a A a u -u + u p+l = b. 

We seek 27r-periodic solutions of (|A-1|) , To this end, first notice that since 
a > ^ one can show by arguments identical that in [FL12, Appendix B] that 
any u G H a (R/2irZ) satisfying (1A7TT) automatically lies in H 2a+1 (R/27rZ). 
Iterating now the identity 

u = [k a {K a + I)}' 1 (u p+1 + {k a - l)u - b) 

implies that u G H°°(R/2ttZ), so that any 27r-periodic H a solution of (lA-lj) 
is automatically a smooth function of z. Thus, it is sufficient to seek solutions 
of (1X11 which lie in i? a (M/2vrZ). 

Set X a := H a (R/2TtZ) x M+ x R, considered to be equipped with the 
natural graph norm, and define the map F : X a — > L 2 (R/2-7rZ) by 

F(v, k, b) := -k a A a v -v + v p+1 - b, 

noting that F is well-defined by Sobolev embedding. Clearly then, the zeros 
of F correspond to 2-7r-periodic solutions of (|A-ip , which can be taken to be 
even functions of z by applying an appropriate spatial translation. First, we 
claim that F is C 1 on X a . Indeed, given any (v , k, b) £ X a we find that 

BF BF 

^— = -k a A a - 1 + (p + l)v p , — = -ak a ^A a 

ov ok 

and ^ = —1. Clearly |jr and ^ depend continuously on (v,k,b) G X a , 
and that ^ depends continuously on (k, b) G R + x R. To see that ^ 
depends continuously on v G if a (R/27rZ), notice that by Sobolev embedding 
we have for all v 1 ,v 2 G H a (R/27rZ) 



n+l p+1 



< C\\V\ - V 2 \\h^ 

= (R/2ttZ) 



I l|P -I- II l|P 

l' Ul ll_ff a (R/27rZ) ' W V2 WH a (B./2irZ) 
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for some constant C > 0, from which we find 

< C\\V\ — V-2,\\lJ<x 



dF dF 

— ( Vl ,k,b)f-—(v 2 ,k,b)f 



X (ll u l|lH«(R/27rZ) + IHItf«(]R/27rZ)) 1 1 / 1 1 H« (R/2ttZ) 

for all / 6 iJ a (IR/27rZ), which establishes the continuity on v, as claimed. 
Together, the above arguments verify that F is C 1 on X a . 

Continuing, by inspection we see that F(Q b ,k,b) = for \b\ <C 1 and 
any k G M + . Furthermore, for fixed 6^1 and k > we have 

— (Q 6 , fc, 6)) = —k a A a - 1 + (p + 1)Q?, 

which, by Fourier analysis, has a trivial kernel in L 2 (]R/27rZ) provided k a / 
k*(b), where fc*(6) is given in the statement of the theorem. When k a ^ 
k* (b) then, the Implicit function theorem implies that the root (Qb,k,b) of 
F(v, k, b) continues uniquely for |A;— ko\ <C 1 and \b— bo\ <C 1. By inspection, 
this must correspond to the nearby equilibrium solutions Qb- To find non- 
constant solutions then, we must consider the case k a = k*(b), in which case 
the kernel of the above linear operator is two dimensional and is spanned by 
e ±tz . We now construct periodic solutions v to the equation F(v,k,b) = 
by using a Lyapunov- Schmidt reduction for \b\ <C 1 and \k a — k*(b)\ <C 1. 
We set k a = k*{b) + k and 

(A-2) v(z) = Q b + X -Ae iz + ^Ae~ lz + h{z) 

where A £ C and h G # a (R/27rZ) satisfies 



/ 



/i(z)e ±i2 dz = 0. 



Substituting these expressions into the equation F(v, k,b) =0 leads to an 
equation of the form 

(A-3) L b h = JV(h,A,A,k,b), 

where 

L b :=^-(Q b ,k*(b),b) 

and A/"(0, 0, 0, fc, 6) = for any k, b G R. Again using Fourier analysis, we see 
that the kernel of L b is two dimensional and is spanned by e ±lz . We denote 
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by P : L 2 (R/2ttZ) — >• ker(Lb) the spectral projection onto the kernel of Lb, 
defined for any u G L 2 (M/2vrZ) by 

Pu(z) = u(l)e iz + u(-l)e~ iz . 

Since Ph = then, the perturbation equation ()A-3p is equivalent to the 
system 



(A-4) 



L b h = (I — P) J\f(h, A, A, k, b) 
= PJV(h,A,A,k,b), 



which is valid for any \b\ <C 1. 

Notice the restriction of Lf, to Z := (Id — P)H a (M./2ir'Z) has a bounded 
inverse given by 



■u(n) 



(1 - |n| Q ) 

for any v £ Z. In particular, this formula shows that the operators (L&l^) 1 
form a family of bounded linear operators depending analytically on b. 
Therefore, the system (|A-4|) is equivalent to 



(A-5) 



v = L- 1 (Id - P) Af(h, A, A, k, c) 
= PAf(h,A,A,k,c). 



Using the implicit function theorem, we can solve the first equation in (|A-5j) 
to find a unique solution h = H*(A, A,k,b) £ Z that depends analytically 
on (A,A,k,b) in a neighborhood of (0,0,0, bo) in the space diag(C 2 ) x M 2 , 
where diag(C 2 ) := {(z,z) : z G C). In particular, notice the uniqueness of 
this solution implies that 

H*(0,0,k,b) = Q b 

for all \b\ <C 1. Furthermore, the invariance of (|2.2p with respect to spatial 
translations z t— > z + zq and the reflection z i— )• — z implies the relations 



(A-6) 



H,(AAM)(^^o) = ^(^e iZ0 ,^e-^°,fc,6) 



Substituting the above into PM(h, A, A, k,c) = now yields the equation 

PAf(iT*(AAM)(2),A^,M = o 
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which must be solved. Using the explicit form of the projection P, this 
equation has solutions provided the two orthogonality conditions 

/•27T a iz i Zp-iz 

J±(A,A,k,b)= / J\f(H*(A,A,k,b)(z),A,A,k,c)dz = 

Jo 2 

are satisfied. Notice that the relations (|A-6|) imply that the functions J± 
satisfy 

J+(Ae te °, Ae"^, k, b) = J + (A, A, k, b) = J + (A, A, k, b) 
J-(Ae iZ0 ,Ae- iZ0 ,k,b) = J_{A,A,k,b) = -J_(A,A,k,b). 

In particular, taking zq = — 2arg(^4) in the equalities for J_ implies that 

J_{A,A,k,b) = J_{A,A,k,b) = -J_(A,A,k,b), 

from which we see that the condition J-(A, A, k, b) = is always satisfied. 

As for the solvability condition J + = 0, taking zq = — arg(A) in (|A-7[) 
we find J+(A, A, k, b) = J+(\A\,\A\,k,b) so that the associated solvability 
condition becomes J+(a, a, k, c) = where a 6 R belongs to a sufficiently 
small neighborhood of the origin. Using the explicit form of jV, we find 

/•27T 

J+(a,a,k,b) = / acos(z)Af(H*(A, A, k, b)(z), A, A, k, c)dz 
Jo 

= a 2 k + J(a, k, b) 

where J is analytic in its argument, even with respect to the parameter a, 
and satisfies J(0, 0, b) = J(0, 0, b) = 0. By the Implicit Function Theorem 
again, for < \a\ 1 we obtain a solution k(a,b) of J + (a,a,k,b) = 
defined for sufficiently small a, b £ R. Furthermore, it follows that k is even 
in a and analytic in a sufficiently small neighborhood of the origin in R 2 , 

so that the function k(a,b) = [k*(b) + k(a,b)j satisfies the properties 
discussed in (i). 

From the above considerations, it follows that the system (|A-5|) has a 
unique solution 

(h,k)=(H*(A,A,k(\A\,b),b),k(\A\,b)), 

defined for any sufficiently small A G C and \b\ <C 1. Substituting h = 
H*(A,A,k(\A\,b),b) into (jA-2j) thus yields a 27r-periodic solution of (f23]l . 
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The periodic solutions P a ^ described in the theorem are now found by re- 
stricting to A G R, so that 

P a ,b(z) = Qb + acos(z) + Va,b(z), v a ,b(z) = ff*(a, a, k(a,b),b)(z) . 

The properties of P a ^ described in (ii) are now easily deduced from ana- 
lyticity and the symmetries of the function H*{A, A, k(\A\, b), 6), while (iii) 
follows directly from the expansion of P a & given in (ii). □ 

B Proof of Proposition 13.11 

In this appendix, we establish a nonlocal type of Floquet-Bloch theory that 
is suitable for the stability analysis presented in this paper. 

Proof of Proposition HOI Clearly (iii) holds for some £ G [—1/2,1/2) if and 
only if the kernel of the operator M a ,b£ — AI acting on L 2 (R/2-7rZ) is non- 
trivial. Furthermore, the operator M. a ,b,£ acting on L 2 (R/27rZ) has densely 
defined and compactly embedded domain ff a+1 (R/2-7rZ) and hence has com- 
pact resolvent. As a result, the spectrum of M a ,b,€ consists of isolated eigen- 
values of finite multiplicity and, in particular, we see that X £ a (M. a ,b,s,) if 
and only if the operator M a ,b,£ ~ AI has a nontrivial kernel. This establishes 
that (ii) <?=^ (iii). 

Now, assume that (ii) does not hold, i.e. that the operator M. a ,b,i ~ AI 
is boundedly invertible on L 2 (R/2vrZ) for all [-1/2, 1/2). Then there exists 
a constant C > such that 

\\(M a ,b£ ~ AI)u|| i 2( R / 27rZ) > C\\v\\ L 2 (r/27tZ) 

for all v G H a+1 (M./2irZ) and £ G [-1/2, 1/2). Using d£TJ then, we find for 
all w G H a+1 (R) that 

rl/2 

\\(M a , b -XI)w\\l 2m =2tt J ^ J(M aM -XI) w(£, -W&w^dt 
>27rC 2 ^ / V(e,-)||i 2(R/27rZ) ^ 

= C 2 |kllL 2 (R)- 

It follows then that M a ,b — AI is boundedly invertible as an operator acting 
on L 2 (R). This establishes that (i)=^(ii). 

Finally, assume that (ii) holds, i.e. that for some £o G [—1/2,1/2) the 
operator A4 a ,b,£ — AI is not boundedly invertible on L 2 (R/2vrZ). Then 
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by the above considerations there exists a v G i? a+1 (]R/27rZ) such that 
{M a ^ - AI) v = 0. For each < e < 1, let (j> £ : M ->■ R be defined by 

e- 1 ' 2 if |f | < e/2 
if \i\ > e/2, 



and note that ||^ e [|l 2 (k) = Given a fixed e S (0, 1) then, notice that the 
function 

M 2 3 &z)^v{z)MZ-Zo)£C 
belongs to L 2 ([-l/2, 1/2); L 2 (IR/2vrZ)), and can hence be viewed as the 
Bloch-transform of some function v £ G L 2 (IR) with 



2 /" 1/2 

IKIIl 2 (K) = J ^ ^eH - ft)) ll w lli2( R /27rZ)^ = ll' 1 ' Hi 2 (M/2ttZ) 

for all e > sufficiently small. Then for each < e <C 1, we find using (|3.1I 
that 

,1/2 

11(^0,6 - AI)« £ ||^ 2(M) = J - Co) 2 IK^W " AI)K-)lli 2(R/2 ^) dC 



- / \\{M aM -\I)v(-)\\ 2 L2W2nZ) dt 



'|«-&l<6/2 

Next, we want to show the above quantity tends to zero as e ->• + . 

To this end, notice that Lemma [3.31 in Section [3.31 implies that the map- 
ping f — > Ai a ,b,^ is continuous in the operator norm from H a+1 (M./2n7i) to 
L 2 (M/2vrZ). Indeed, for a given w G F Q+1 (M/27rZ) and fi,6 G [-1/ 2 , 1/2) 
we have the estimate 

\\(M a ,b,Si ~ ■ M a,b,&)w\\ L 2( R/27T Z) ~ ICl - 62 1 1 1 ^ | [ ifa+i (R/2ttZ) j 

yielding the desired continuity. It now follows that 

7 / \\(M aM - AI) uCOIliaon/M) d ^ 

= \\(M aMo - AlHli2 (K/27rZ) = 0, 
so that, in particular, given any n»l there exists an e n > such that 

2 1 

\\{M a , b - XI)v £ J L2m < -. 

Recalling that ||w £ ||l2( K ) = |M|l 2 (k/27tZ) for all < £ <C 1, it follows that the 
operator M. a ,b — AI is not boundedly invertible. This establishes (ii)=^(i), 
which completes the proof. □ 
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